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An  evaluation  of  a  series  expansion  method  for  predicting  the  far- 
field  from  near-field  cross-power  measurements  of  partially  coherent 
sound  fields  is  given.  This  series  method  is  based  on  an  expansion  in 
spherical  wave  functions  truncated  to  include  only  the  number  of  terms 
equal  to  the  number  of  near  field  measuring  points  (which  need  not  be 
restricted  to  the  usual  separable  surfaces).  The  evaluation  is  through 
a  set  of  numerical  examples  where  the  predicted  far  field  is  compared 
with  the  exact  far  field  from  a  theoretical  source  of  partially 
coherent  sound.  It  is  shown  that  good  results  can  be  obtained  for 
incoherent,  partially  coherent,  and  coherent  sound  fields. 
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I.  Introduction 

■  la  this  reporC  we  shell  demonstrate,  through  e  numerical  evalua¬ 
tion,  that  the  far  field  can  be  predicted  from  the  near  field  of 
coherent,  partially  coherent,  and  incoherent  sources  by  a  finite  series 
method.  The  series  method  has  been  described  in  some  detail  in  an 
earlier  memorandum^  where  an  evaluation  of  the  method  for  only  coherent 
sources  was  given.  In  order  to  test  the  method  for  various  degrees  of 
coherency  it  was  necessary  to  devise  a  theoretical  model  for  a  parti¬ 
ally  coherent  source  distribution.  The  details  of  this  model  are  the 

2 

subject  of  another  technical  memorandum  . 

The  numerical  evaluation  presented  in  this  memorandum  is  based  on 
a  comparison  of  the  exact  and  predicted  far  field  intensity  distribu¬ 
tion.  The  exact  far  field  is  calculated  directly  from  the  theoretical 
source  distribution  while  the  predicted  far  field  Is  obtained  by  the 
series  method  using  near  field  calculations  of  the  cross  power  from  the 
theoretical  source  distribution.  These  near  field  samplings  simulate 
the  near  field  measurements  that  would  be  made  by  the  hydrophones  about 
an  actual  source. 

Ve  have  confined  our  attention  in  this  report  to  distributions 
which  have  axial  synutetry.  A  brief  review  of  the  series  method  is 
given  along  with  a  brief  description  of  the  theoretical  model  of  the 
partially  coherent  source  distribution  used  to  test  the  method.  The 
results  of  the  far  field  predictions  are  then  compared  with  the  exact 
values  for  various  hydrophone  and  source  distributions. 
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[I.  Review  of  the  Series  Method 


If  we  consider  axially  symmetric  sources,  Che  solution  for  the 
cross-power  function  may  be  written  as 

C  2  -Z  ($r£)  rh(c**i )  ltm  (£>y)  £  (cac6j  )  M 

where  the  conventional  notation  for  spherical  wave  functions  has  been 
adopted,  and  the  are  the  coefficients  to  be  determined.  It  is 

assumed  that  the  fields  of  Interest  may  be  adequately  represented  by  a 
finite  number  of  terms,  H  .  Consider  now  ti  measuring  hydrophones 
(which  need  not  be  on  the  usual  separable  surfaces).  The  hydrophones 
would  measure  the  auto  spectral  density  (intensity)  V  .  ,  w) 

at  the  H  positions, and  the  cross-power  u> )  between  all  the 

pairs  of  the  N  hydrophones  can  be  determined  yielding  a  total  of 
measurements  at  each  frequency.  With  the  two  summations  in  the  above 
equation  each  truncated  to  Include  only  the  first  terms  there  will 
then  be  /V4,  coefficients  which  can  be  determined. 

Under  the  conditions  of  only  /V  terms  in  the  sums,  the  above 
equation  may  be  written  more  compactly  in  the  matrix  notation. 

Adopting  the  summation  convention  the  above  may  be  written  (in  terms  of 
the  elements)  as 

-  K.  a.mi  <2> 


where 


4,;-  C 
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and  the  asterisk  means  complex  conjugate.  With  the  upper  case  letters 
denoting  the  corresponding  matrices  we  have.  In  matrix  notation 

«-A*CA  (4) 

where  A*  Is  the  transpose  of  the  complex  conjugate  of  A  .  The  solu¬ 
tion  for  the  coefficients  can  then  be  obtained  through 

matrix  Inversion  so  that 

C-  A*  G  A'*  .  (5) 

# 

A  solution  for  the  entire  field  Is  obtained  with  the  evaluated 
substituted  Into  the  expansion.  The  degree  to  which  the  result  Is  an 
accurate  prediction  depends  on  the  number. and  location  of  the  hydro¬ 
phones. 

III.  Description  of  the  Source  Model 

The  development  of  the  source  model  Is  given  In  detail  In  an 
2 

earlier  memorandum  .  The  model  specialized  for  axial  symmetry  consists 
of  a  linear  distribution  of  small  spherical  acoustical  radiators;  that 
is,  point  sources.  The  time  dependence  for  the  nC^  source  Is  assumed 
to  be  a  general  function  •  Tht  solution  for  the  pressure  fa  at  fj 

can  then  be  written  as 

.  g  WO 

<6> 

where  M  is  the  number  of  sources,  V*.  is  the  distance  from  the 
source  to  the  field  point  at  **•  ,  and  C  is  the  velocity  of  sound. 

The  cross-power  S,^C ^ ,  r .  ,  u>)  in  the  field  was  obtained  from  the 
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product  of  the  Fourier  transform  of  ^  at  vj  ,  end  the  complex  conju¬ 
gate  of  the  transform  at  r.  .  On  aeparatlng  the  factors  which  depend 
on  the  sources  only  and  those  which  relate  the  sources  to  the  field  the 


result  may  be  written  in  the  form 


m  M 


n-  ri«  <7) 

where  is  the  cross  power  amongst  the  sources.  By  writing  the 

expansion  in  this  way  we  need  not  specify  the  exact  phase  of  the 
sources  and,  in  fact,  they  may  be  random.  What  is  specified  is  the 
intensity  of  each  source  and  the  cross-power  or  correlations  between 
the  sources.  Thus,  for  a  chosen  (**>)  the  near  field  as  well  as  the 
far  field  cross  £|t  may  be  calculated  through  Eq.  (7)  for  Yi  sources 
at  frequency  tO . 

We  would  like  to  emphasize  that  this  calculation  does  not  imply 
the  sources  are  necessarily  single  frequency  (and  thus,  coherent). 

What  is  meant  by  a  calculation  at  a  single  frequency  is  that  this  is 
one  frequency,  of  a  possible  multitude,  at  which  we  focus  our  attention 
In  actual  practice  this  can  be  done  by  frequency  filters  which  ideally 
do  not  change  the  degree  of  coherency. 

One  of  the  difficult  problems  is  to  choose  a  representative  q. 

As  an  initial  test  we  chose 

-  ) 

.  *  /  A  \ 


where  we  assumed  the  sources  to  be  spaced  an  equal  distance  5  .  The 
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quantity  clw-nl  represents  the  separation  between  the  sources  and 
we  call  j£  the  coherence  length.  For  ^  »  S  I  m-n  (  it  can  be  seen 
that  is  unity  and  all  the  sources  are  coherent,  while  for 

£  «.  S  !.n-w\  vanishes  for  v»^L>w  and  the  sources  are 

incoherent.  For  intermediate  values  of  £.  the  sources  are  partially 
coherent. 

IV.  Numerical  Evaluation 

Here  we  present  some  of  the  results  of  a  numerical  evaluation  of 
the  series  method.  The  calculations  were  programmed  for  computation  on 
our  IBM  1130.  Because  of  limited  core  storage  it  was  convenient  to 
evaluate  arrays  corresponding  to  only  eleven  hydrophones.  The  program 
for  the  numerical  evaluation  can  be  separated  into  three  parts: 

1)  The  "source  program"  -  from  which  the  exact  near  and  far  field 
values  are  calculated  using  various  source  distributions  described  by 
equations  (7)  and  (8). 

2)  The  "coefficient  program"  -  from  which  the  G  coefficients  are 
evaluated  using  matrix  equation  (5)  along  with  the  near  field  cross 
power  values  from  the  "source  program". 

3)  The  "field  prediction  program"  -  from  which  the  far  field  is 
calculated  using  equation  (l)  along  with  the  coefficient  values  obtained 
from  the  "coefficient  program". 

The  scheme  for  the  display  of  the  results  may  be  understood  by 
referring  to  Fig.  1.  The  separation  and  relative  location  of  the 
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sources  and  bydrcphones  (points  at  which  the  near-field  values  are  used) 
is  shown  piccorially  in  the  upper  left-hand  quadrant.  Aa  shown  here, 
the  evaluation  is  for  three  sources  each  separated  by  one-half  wavelength 
(X)  with  near-field  values  sampled  at  the  eleven  hydrophone  locations. 
The  hydrophone  array  is  symmetrical  and  at  a  normal  distance  of  1.59 
wavelengths  with  inter-element  spacing  given  by  0.318  wavelengths. 

The  results  for  three  degrees  of  coherency  are  shown  in  the 
remaining  quadrants  of  the  page.  The  exact  values  for  the  far  field 
intensity  £|(  £?-,  t‘  ,  are  given  as  a  solid  line  and  the  predicted 
values  are  given  by  the  dots.  The  three  degrees  of  coherency  were 
obtained  using  Eq.  (8)  with  the  values  of  Jljx  equal  to  0.1,  1.0,  and 
10.0;  and  calling  these  cases  incoherent,  partially  coherent,  and 
coherent,  respectively.  Except  for  Fig.  2  coherency  is  essentially 
reached  for  /.Is.  equal  to  ten. 

The  results  for  the  sources  have  been  normalized  so  that  each 
source  separately  yields  an  amplitude  value  unity  in  the  far  field. 

Thus,  for  three  incoherent  sources  we  have  an  omnidirect? cnal  pattern 
with  an  intensity  value  of  three  since  here  the  intensity  is  simply  the 
sum  of  the  intensities  of  each  source.  For  complete  coherency  the  far 
field  intensity  result  at  90°  ia  Che  square  of  the  sum  of  amplitude 
values  yielding  a  value  9  for  three  sources;  moreover,  we  see  that 
here  we  get  the  usual  pattern  structure  with  side,  lobes  and  nulls.  The 
case  of  partial  coherency  is  characterized  by  a  main  lobe  without  deep 
nulls  and  separate  side  lobes.  The  exact  values  for  the  far-field 
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cross-power  for  this  source  distribution  ere  given  in  en  earlier  ceio* 

2  o 

rand urn  (in  which  the  coordinate  system  was  rotated  90  relative  to  the 

one  in  this  memorandum). 

Inspection  of  the  predicted  and  exact  far-fleld  intensity  results 
shown  in  Fig.  1  indicate  quite  clearly  that  for  this  case  good  agree* 
ment  is  obtained  for  all  three  degrees  of  coherency.  In  fact,  for  the 
case  of  Incoherency  the  predicted  values  of  the  intensity  at  90° ,  45° 
and  0°  are  3 *00001,  3*00028  and  3*00249  respectively.  If  we  compare 
the  worst  predicted  value  at  0°  (l.e.,  along  the  axis  of  the  source) 
with  the  exact  value  3j  Che  worst  percentage  error  is  less  than  one- 
tenth  of  one-percent.  This  very  good  agreement  is  partly  due  to  the 
large  number  of  hydrophones  (l.e.,  near  field  sampling  points)  compared 
to  the  ratio  of  the  overall  source  length  to  wavelength.  For  best 
results  we  would  expect  (for  reasonable  distributions  and  array  loca¬ 
tions)  that  the  number  of  hydrophones  exceed  5(  wX  )  where  L  is  the 
length  of  the  hydrophone  line.  In  this  particular  case  L  is  equal  to 
X  and  the  eleven  hydrophones  safely  exceed  the  number  five. 

In  Fig.  2,  the  hydrophone  arrangement,  overall  source  length,  and 
£/s  ratios  are  the  same;  however,  here  eight  additional  sources  have 
been  Inserted  between  the  two  outer  sources.  In  this  case  the  separa¬ 
tion  between  sources  is  one-tenth  of  one-wavelength,  and  the  source 
acts  as  a  nearly  continuous  line  one-wavelength  in  length.  Again  the 
agreement  between  the  exact  and  the  predicted  values  of  the  intensity  is 
good,  and  the  added  complication  of  more  sources  does  not  appear  to 
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product  any  aarioua  problems  providing  that  tha  ovarall  length  of  the 
sourca  line  la  not  Increased,  tfa  should  note  that  McoharencyN  has  not 
baen  actually  reached  in  this  casa  sine «  the  coherence  length  is  now  on 
the  order  of  the  overall  source  length  rather  than  such  greater  than 
the  sourca  length  as  it  wss  in  Fig.  1. 

In  Fig.  3  sane  hydrophone  array  is  used  at  the  sane  location; 
however,  hare  three  sources  are  used  again,  and  this  tine  the  spacing 
is  one  wavelength  for  an  overall  length  of  two  wavelengths.  Va  see 
that  in  all  three  cases  tha  prediction  is  quite  good  in  the  arc  fron 
50°  to  90°.  This  night  be  expected  since  this  region  is  covered  quits 
adequately  by  the  hydrophone  array.  The  results  are  for  the  nost  part 
in  poor  agreement  in  tha  remaining  region  for  the  coherent  case; 
however,  we  see  that  there  is  a  trend  toward  better  agreement  in  the 
axial  direction  (0°)  as  the  coherency  decreases. 

For  the  case  under  discussion  the  number  of  hydrophones  (11)  is 
approximately  equal  to  rather  than  considerably  greater  as  in 

the  previous  cases.  The  results  can,  however,  be  improved  by  increas¬ 
ing  the  separation  between  the  hydrophones.  This  improvement  is  shown 
in  Fig.  4  where  now  the  separation  between  hydrophones  is  one  wave¬ 
length  whereas  in  the  previous  cases  it  was  approximately  one-third  of 
',ue-wavelength.  The  prediction  deteriorates  as  the  hydrophone 
separation  increases  much  beyond  one  and  one-quarter  wavelengths. 

Tha  results  for  a  shorter  hydrophone  array  in  closer  to  the  three 
sources,  with  source  separation  X/l  ,  are  shown  in  Figs.  5  and  6.  The 
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separation  between  the  hydrophone  elements  is  new  .064  X  so  that  the 
total  length  of  the  hydrophone  array  is  less  than  the  overall  length  of 
the  source.  In  Fig.  5  the  normal  distance  to  the  hydrophone  array  is 
slightly  beyond  the  smallest  circle  drawn  through  the  two  outer  sources 
whereas,  in  Fig.  6  the  same  array  is  close  enough  so  that  all  the 
hydrophones  are  entirely  within  this  circle. 

tfe  see  that  the  hydrophone  line  is  too  close  and  too  short  for  an 
accurate  prediction  throughout  the  field.  Although  the  more  remote 
line  location  gives  better  results  in  the  direction  of  the  main  lobe, 
the  line  location  within  the  circle  yields  a  better  average  prediction. 
It  appears,  therefore,  that  a  reasonable  prediction  of  the  far  field 
can  be  obtained  even  if  all  the  hydrophones  are  within  the  smallest 
circle  through  the  outer  sources.  Satisfactory  agreement  has  also  been 
obtained  for  cases  where  part  of  the  hydrophone  arr*y  was  within,  and 
the  remaining  part  outside  of  this  circle.  The  fact  that  a  reasonable 
prediction  can  be  obtained  even  though  the  sampling  points  (hydrophone 
locations)  are  within  this  circle,  is  of  theoretical  as  well  as  practi¬ 
cal  interest  and  requires  further  study. 

V.  Summary  and  Conclusions 

The  series  expansion  method  for  the  prediction  of  the  far  field 
from  near  field  measurements  of  coherent,  partially  coherent,  and 
Incoherent  acoustic  fields  has  been  evaluated  through  a  set  of 
numerical  examples.  It  was  found  that  good  agreement  can  be  obtained 
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for  *11  degrees  of  coherency.  The  shortcomings  of  this  series  method 
for  cross-power  measurements  sppesr  to  be  no  worse  than  those  found  for 
the  series  method  for  pressure  amplitude  and  phase  measurements  reported 
earlier^. 
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